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Abstract 



A non-LTE analysis of K I resonance lines at 7664.91 and 7698.97 A was carried out for 15 red giants 
belonging to three globular clusters of different metallicity (M 4, M 13, and M 15) along with two reference 
early-K giants (p Boo and a Boo), in order to check whether the K abundances are uniform within a cluster 
and to investigate the behavior of [K/Fe] ratio at the relevant metallicity range of —2.5 < [Fe/H] < — 1. 
We confirmed that [K/H] (as well as [Fe/H]) is almost homogeneous within each cluster to a precision of 
| ^0.1 dcx, though dubiously large deviations are exceptionally seen for two peculiar stars showing signs 

O |. of considerably increased turbulence in the upper atmosphere. The resulting [K/Fe] ratios are mildly 

supersolar by a few tenths of dex for three clusters, tending to gradually increase from ~ +0.1-0.2 at 
[Fe/H] ~ -1 to ~ +0.3 at [Fe/H] - -2.5. This result connects reasonably well with the [K/Fe] trend 
of disk stars (—1 < [Fe/H]) and that of extremely metal-poor stars (—4 < [Fe/H] < —2.5). That is, 
[K/Fe] appears to continue a gradual increase from [Fe/H] ~ toward a lower metallicity regime down 
to [Fc/H] ~ —3, where a broad maximum of [K/Fc] ~ +0.3-0.4 is attained, possibly followed by a slight 
downturn at [Fc/H] < -3. 
■ Key words: stars: abundances (potassium) — stars: atmospheres — 

stars: globular clusters (M 4, M 13, M 15) — stars: late-type — stars: spectra 
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■<^j- ■ 1. Introduction lar atmospheric parameters as well as the line strength. 



Actually, since the extent of the (negative) non-LTE cor- 



The chemical evolution of potassium (K) in the Galaxy, rection amounts from a few tenths dex even up to ~ 1 dex, 

' which can be investigated by examining its photospheric it is requisite to take account of the non-LTE effect in de- 

. abundances of old metal-poor stars, is still only insufh- riving stellar K abundances. 

ciently explored and not yet well understood, especially in — Besides, they occasionally suffer serious blending with 

the metal-poor regime of halo stars. The main reason for strong absorption lines originated from earth's atmo- 

1 this unsatisfactory situation is presumably the difficulty sphere (especially for the 7665 line) because of being lo- 

\ involved with determinations of stellar potassium abun- catcd in an unfavorable wavelength region crowded with 

dances, for which only two strong K I lines at 7664.91 and such telluric lines. Hence, when encountered with such 

7698.97 A (hereinafter referred to as 7665 and 7699 lines unfortunate cases, one has to recover the pure stellar spec- 

for brevity) are practically usable. That is, although these trum by appropriately eliminating the blended telluric 

resonance lines are so strong as to be measurable even for spectrum. 

stars of extremely low metallicity, the following features The first non-LTE study of K abundances in metal- 
make the problem comparatively harder to deal with: poor stars was carried out by Takeda et al. (2002a, here- 
- Above all things, they generally suffer an apprccia- 

ii i rmij a j. 1 i*i i l ,i ,i carried out a detailed study on the formation of this line in 

bly large non-LlE effect, which depends on the stel- „ , , ... „ + , , „, . , .„ , n „„„ N 

J ° r Procyon (along with the bun), ivanova and bnimanskn (2000) 

* Based on data collected at Subaru Telescope, which is operated calculated non-LTE abundance corrections applicable to A-K 

by the National Astronomical Observatory of Japan. stars of wide Parameter ranges. Similarly, Takeda et al. (2002a) 

+ rp. i . ■ i ui -ci -ii u j -l ui i 4-u r>Ao7 carried out extensive non-LTE calculations on a grid of models 

1 the electronic table El will be made available at the PASJ a 

i ., uv ,. i ■ i • i ■ ■ ,, , , . and published tables of non-LTE corrections applicable to F— G— 

web site upon publication, while it is provisionally placed at 1 , 

(http://optik2.mtk.nao.ac.jp/-takeda/GCpotassium/). K dwarfs throu S h ™pergi a nt S . Recently, Zhang et al. (2006a) 

1 While non-LTE line formation of K I 7665/7699 resonance lines performed a careful non-LTE investigation of solar potassium 

■ n.jv i i u uji, ij-i n«n i lines with a special attention to clarifying the important atomic 

specific to the solar atmosphere had begun already in 1980s (see 1 J a f 

-a , . , 1nnn j 4.1. c j-u • 1 x rjiT-i . parameter of neutral hydrogen collision cross section, which they 

Bruls et al. 1992 and the rctercnccs therein), non-LlE studies *; jo > j 

r .1 -*/■ T a I, . v .■ .I. , ,i ,i ,i o further applied to determinations of K abundances for mctal- 

ot these K I doublet lines directed to stars other than the Sun *+ 

gradually appeared in the last decade: Takeda et al. (1996) poor stars ( Zhan S & Zhao 2005 ' Zhan S et aL 2006b )- 
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inafter referred to as Paper I). They showed that [K/Fc] 
values for the main-sequence stars in the galactic disk 
(— 1 < [Fe/H] < 0) tend to show a rather tight relation 
of steadily increasing with a decrease of metallicity (such 
as like the trend of a-group elements). This was actu- 
ally a reconfirmation of the results of Chen et al.'s (2000) 
LTE analysis, despite the significant non-LTE corrections 
of ~ —0.5 dex, because the corrections turned out to act 
similarly on the Sun and these disk dwarfs (i.e., [K/H], 
the differential stellar potassium abundance relative to the 
Sun. was eventually not much affected). Such an a-like 
trend of [K/Fe] for disk stars was also confirmed by Zhang 
et al. (2006b). 

Unfortunately, the situation becomes quite uncertain 
when we go into the much lower metallicity regime 
([Fe/H] < — 1). It was already noticed in Paper I from 
the non-LTE reanalysis of the published data for several 
halo stars of [Fe/H] - -2 (Gratton & Sneden 1987b) that 
their [K/Fe] values widely spread from ~ +0.1 to ~ +0.7 
(cf. figure 4a of Paper I; see also figure 8 of Gratton & 
Sncdcn 1987a). Zhang and Zhao (2005) also suggested a 
considerably large diversity of [K/Fc] (from ~ to ~ 1) 
for the halo stars of — 2 < [Fc/H] < — 1 they studied (cf. 
figure 12 therein). On the other hand, however, five halo 
stars included in Zhang et el.'s (2006) study have nearly 
the same [K/Fe] ratios around ~ +0.2 over the range of 
— 2< [Fc/H] J$ — 1 (cf. figure 5 therein). Besides, Cayrelet 
al. (2004) obtained in their extensive study of extremely 
metal-poor stars (-4< [Fc/H] < -2.5) that [K/Fc] shows 
a slightly supersolar trend (~ +0.1-0.2) decreasing toward 
a very low metallicity with a fairly small scatter (cf. figure 
9 therein). Thus, at present, we know little about how the 
[K/Fe] ratio actually behaves itself at [Fe/H] < -1: Tight 
trend? Large diversity? Does it smoothly connect to the 
a-like clear tendency seen in disk stars at [Fe/H] > -1? 

For the purpose of checking the possibility of K abun- 
dance diversity in such metal-deficient halo stars, it would 
be interesting to study stars in globular clusters, for which 
any trial of potassium abundance determination has never 
been reported so far to our knowledge. Admittedly, the 
elemental abundances of old globular cluster stars may 
not necessarily be the same as those of field halo stars, 
given that specific abundance peculiarities are known for 
several elements (e.g., O, Na, Mg, Al) which are presum- 
ably due to dredge-up of nuclear-processed product caused 
by evolution-induced deep mixing. However, since K is 
considered to be synthesized mainly via the oxygen burn- 
ing in high-mass stars, whichever the process is explosive 
or hydrostatic (see, e.g., tabic 19 in Woosley & Weaver 
1995), it is unlikely that the surface K abundance under- 
goes any appreciable changes during the course of stellar 
evolution in low-mass stars of globular clusters. Besides, 
the fact that no star-to-star variation in globular clusters 
is observed in the abundances of Ca [Z = 20, near to 
Z = 19 of K) may also suggest the inertness of K to nu- 
clear processes in the stellar interior, presumably because 
of the higher Coulomb barrier compared to lighter ele- 
ments. Therefore, as a reasonable working hypothesis, we 
may postulate (such as the case for Fe in most clusters) 



that (i) the abundance of K was the same in any stars 
of a given globular cluster when they were born, that (ii) 
this uniformity has been retained up to now, and that (iii) 
the surface K abundances of cluster stars may be regarded 
to be equivalent to field halo stars of similar metallicity. 
Accordingly, if we could observationally confirm postula- 
tion (ii), this would assure the practical validity of (iii), 
which means that we may directly compare the results of 
cluster stars with those of other halo stars in general. Or, 
alternatively, if we found a markedly large diversity of K 
abundances within a cluster contrary to postulation (ii), 
we would have to consider a possibility of real star-to-star 
variation (or cast doubt on the validity of our abundance 
determination method) . 

Motivated by this consideration, we decided to conduct 
a spectroscopic study on the intrinsically bright red-giant 
stars of three globular clusters (M 4, M 13, and M 15) 
covering a wide metallicity span (from [Fe/H] ~ —2.5 to 
~ — 1) based on the high-dispersion spectra obtained with 
Subaru/HDS, in order to determine the abundance of K 
for each star while taking into account the non-LTE effect, 
after having established the atmospheric parameters in- 
cluding [Fe/H]. What we intend to clarify is, as described 
above, to check the homogeneity of [K/H] within a clus- 
ter [postulation (ii)], and to study the behavior of [K/Fc] 
in halo stars based on the results of these globular clus- 
ters [postulation (iii)] while comparing with those of other 
metallicity regime. This is the purpose of this study. 

Besides, in connection with the main subject, we car- 
ried out a reanalysis of the equivalent-width data of K I 
7665 and 7669 lines for ~ 30 extremely metal-poor stars 
published by Cayrel et al. (2004), in order to establish 
the [K/Fe] vs. [Fe/H] relation at the metallicity range of 
—4 < [Fe/H] < —2.5. The main motivation is to properly 
take into account the non-LTE effect, since their treat- 
ment in this respect does not appear to be sufficiently 
valid (i.e., they applied a tentatively assumed correction 
uniformly to all stars). This reanalysis is separately de- 
scribed in the Appendix. 

2. Observational Data 

Three globular clusters (M 4, M 13, and 
M 15) were chosen for this study, because (1) 
they have different metallicities from each other 
([Fe/H] M i 5 < [Fe/H] M13 < [Fe/H] M4 ), (2) they are 
located in near or reasonable distance to us (so that red 
giants of y ^ 12 mag exist), and (3) they are compara- 
tively well studied and a number of literature data may be 
found. Practically, we selected 5 red-giant stars for each 
cluster satisfying the criterion of 4100 K < < 4300 K 
(j>p ° i s t ne effective temperature photometrically evalu- 
ated from (B — V)o colors; cf. subsection 5.1): L 2406, 
L 2617, L 3209, L 3624, and L 4511 for M 4; 1-13, 11-76, 
111-52, 111-59, and 111-73 for M 13; K 144, K 341, K 431, 
K 634, and K 825 for M 15. These targets belong to the 
brightest-class group in each cluster (V ~ 11-12 mag for 
M4,y~ 12-13 mag for M 13 and M 15). In addition, two 



No. ] 

early-K giant stars of near-solar or subsolar metallicity 
(p Boo and a Boo) were also included as the standard 
stars, though they are normal K giants of luminosity class 
III with log(L / Lq) ~ 2 while all 15 cluster stars are "tip 
giants" locating on the tip of asymptotic giant branch at 
log(L/L ) ~ 3. 

The observations of these 17 stars (along with a rapid 
rotator Altair as a reference of telluric lines) were car- 
ried out on the night of 2008 August 20 (Hawaii Standard 
Time) by using the High Dispersion Spectrograph (HDS; 
Noguchi et al. 2002) placed at the Nasmyth platform 
of the 8.2-m Subaru Telescope, which can record high- 
dispersion spectra covering a wavelength portion of ~ 
2600A (in the red cross disperser mode) with two CCDs 
of 2Kx4K pixels at a time. With the slit width set at 0."6 
(300 pm) and a binning of 2x2 pixels, the resolving power 
of the obtained spectra is i?~ 60000. In the standard "Ra" 
setting, our spectra cover the wavelength region of 5100- 
6400 A (blue CCD) and 6500-7800 A (red CCD), which 
was so chosen as to make use of the yellow-red region 
of A > 5000A (where the sensitivity of CCD is large and 
many Fe lines usable for parameter determinations exist) 
while including the targeted K I 7665/7699 lines amply. 
The seeing size was 0."5-0."6, and all of the targets could 
be successfully observed without any significant influence 
of neighboring stars. The integrated exposure time for 
each cluster star was from 10 min (= 5 min x2; for M 4 
stars) to 15-20 min (= 7.5-10 min x2; for M 13 or M 15 
stars). 

The reduction of the spectra (bias subtraction, flat- 
fielding, scattered-light subtraction, spectrum extraction, 
wavelength calibration, continuum normalization) was 
performed by using the echelle package of the software 
IRAF 2 in a standard manner. The resulting average S/N 
ratios were around ~ 100-150 for each of the 15 clus- 
ter stars, while much higher values were attained for the 
bright reference stars; ~ 300-400 (p Boo), - 200-300 
(a Boo), and - 600-700 (Altair). 

3. Atmospheric Parameters 

The determination of the atmospheric parameters [T e ff 
(effective temperature), logg (surface gravity), Vt (mi- 
croturbulence), and [Fe/H] (= Af™ - A%; differential Fe 
abundance relative to the Sun, where Ap c is 7.50 in the 
usual normalization of A-n = 12.00)] necessary for con- 
structing model atmospheres was implemented by way of 
the spectroscopic approach using the equivalent widths 
(EWs) of Fe I and Fe n lines, which has a merit of es- 
tablishing these four parameters based only on the same 
spectrum to be further used for abundance determina- 
tions. 

Practically, we used the computer program (named 
TGVIT) developed for this purpose (Takeda et al. 2005; 
cf. section 2 therein), which is based on the principle 

2 IRAF is distributed by the National Optical Astronomy 
Observatories, which is operated by the Association of 
Universities for Research in Astronomy, Inc. under cooperative 
agreement with the National Science Foundation. 
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of searching for the most optimum solution in the 3- 
dimcnsional (T c g , log<?, v t ) space such that simultaneously 
satisfying the three requirements: (1) %i ow -independence 
of Fe I abundances (excitation equilibrium, where xiow is 
the lower excitation potential), (2) equality of the mean 
abundance derived from Fe I lines and that from Fe II 
lines (ionization equilibrium), and (3) EW-independence 
of the abundances (curve-of-growth matching), as de- 
scribed in Takeda, Ohkubo, and Sadakane (2002b). Since 
low-gravity K-type giants are specifically involved in this 
study, we newly computed a grid of data files [covering 
the parameter ranges of 3750-5000 K in T oS (K), 0.0-2.5 
in \ogg (cm s~ 2 ), 0.0-4.0 in v t (km s _1 ), and -2.5 to +0.1 
in [Fe/H] (dex)] to be used in application of TGVIT by 
interpolation (or extrapolation). 

First, we measured the EWs of available lines by con- 
sulting the list of 330 Fc lines (cf. electronic table El 
in Takeda et al. 2005) by using the Gaussian fitting 
method, where only lines weaker than 100 mA were used 
in order to make sure that errors caused by damping 
wings/parameters are suppressed to a negligible level. 
Then, given these EW data as inputs, the converged so- 
lutions of T ff, logg, and vt (along with Ap e as a product) 
were established by iteratively running TGVIT (cf. sub- 
section 3.2 in Takeda et al. 2005). The resulting parame- 
ter solutions are summarized in table 1, while the detailed 
EW data and the Fe abundances corresponding to the fi- 
nal parameters for each star are presented in electronic 
table El. The trend of the Fc abundances corresponding 
to the final solutions is plotted against EW as well as %iow 
in figure 1, where we can see that the required conditions 
are reasonably accomplished. 

The internal statistical errors involved with these so- 
lutions of T e ff, \ogg, vt, and [Fc/H], which were derived 
by the procedure described in subsection 5.2 of Takeda, 
Ohkubo, and Sadakane (2002b), turned out to be ~ 10- 
30 K, ~ 0.05-0.1 dex, - 0.1-0.4 km s" 1 , and - 0.03- 
0.1 dex, respectively. We will discuss in subsection 4.1 
these spectroscopic parameters in comparison with those 
derived from the conventional method (i.e., use of photo- 
metric colors or evolutionary tracks). 

4. Analysis of Potassium Lines 

4-1- EW Measurement 

Among our 17 target stars, 5 stars of M 13 were found 
to be associated with the most unfortunate case. The K I 
7699 line could not be measured at all because it hap- 
pened to fall just on the narrow inter-order gap of HDS 
(at 7687-7694 A) due to M 13's large (negative) radial 

velocity of 220 to -230 km s~\ while the K I 7665 

line was seriously blended with the strong telluric O2 line 
at 7659.3 A. Such a contamination in the K 1 7665 re- 
gion is also seen for the case of two standard stars, a Boo 
and (especially) p Boo, which turned out to be influenced 
by the telluric 2 doublet at ~ 7664/7665 A. Wc, there- 
fore, tried to eliminate the effect of telluric lines for these 
7 stars by dividing the raw spectra by the spectrum of 
Altair (rapid rotator) by using the IRAF task telluric. 
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Fig. 1. Fc abundance vs. equivalent width relation (left 
panel) and Fe abundance vs. lower excitation potential re- 
lation (right panel) corresponding to the finally established 
atmospheric parameters of T e fj, logg, and vt for each of the 
17 stars. The filled (blue) and open (red) circles correspond 
to Fe I and Fe II lines, respectively. The results for each star 
are shown relative to the mean abundance indicated by the 
horizontal dashed line, and vertically offset by 1.0 relative to 
the adjacent ones. 

The resulting as well as the original spectra are shown 
in figure 2, where we can see that this elimination surely 
worked (though not perfectly). Figure 3 displays the spec- 
tra used for measurements of EWre65 and -EVF7699 , which 
were measured by the Gaussian fitting method as for the 
case of Fe lines. The finally adopted values of EW766S 
and EW7G99 are summarized in table 1. 



Actually, we had to apply further corrections to EW^eas for five 
M 13 stars (1-13, 11-76, 111-52, 111-59, and III- 73) and p Boo, for 
which we tried elimination of telluric lines (EWjae^ for a Boo 
does not need to be corrected, since the important line core was 
free from contamination). Namely, given that this removal could 
not be completely done, as can be recognized from the appre- 
ciable residual at the position of the redward component of the 
relevant 2 doublet (at 7660.4 A for M 13 stars and at 7665.9 A 
for p Boo; see the depression marked with "x" in figure 2), we 
evaluated its equivalent width ( AEW Ies id) and subtracted from 
the directly measured EW as EWlfij^ = EW^ 5 - AEW Ieaid 
for these six stars, assuming that the residuals of both telluric 
components are almost the same. The applied AEW Ics id cor- 
rections arc given in the caption of table 1 . 
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Fig. 2. Removal of telluric O2 lines in the K I 7665 line 
region for the seven stars significantly suffering this effect. 
Dividing the raw stellar spectrum (middle, black) by the spec- 
trum of a rapid rotator Altair (bottom, red) results in the 
final spectrum (upper, blue) where the effect of the blending 
component is reasonably (even if not completely) eliminated. 
The residual caused by the imperfect elimination for the red- 
ward component of the telluric doublet lines is marked by a 
cross (x), the strength of which is usable for the correction of 
EW(K I 7665) similarly affected by the incomplete removal 
of the blueward compenent. The position of the K I 7665 line 
is indicated by an arrow. Each spectrum (normalized with 
respect to the continuum) is vertically offset by 1.0 relative 
to the adjacent one. No Doppler correction has been applied 
to the wavelength scale. 



4-2. Abundance Determination 

The determination of potassium abundances from the 
measured EWrees and EW^gg by taking into account the 
non-LTE effect was carried out in almost the same man- 
ner as in Paper I. As the basic grid of model atmospheres, 
we used Kurucz's (1993) ATLAS9 models (corresponding 
to vt = 2 km s _1 ), which were three-dimensionally inter- 
polated (or extrapolated in special cases of negative logg) 
with respect to T e g, \ogg, and [Fc/H] to generate the at- 
mospheric model for each star. 

The non-LTE statistical-equilibrium calculations were 
implemented for a grid of 90 models resulting from com- 
binations of three T cS values (4000, 4250, 4500 K), five 
logs values (0.0, 0.5, 1.0, 1.5, 2.0), and six [Fe/H] values 
(-2.5, -2.0, -1.5, -1.0, -0.5, 0.0), so that we can obtain 
the depth-dependent non-LTE departure coefficients for 
any star by interpolating (or extrapolating) this grid. See 
Takeda et al. (1996) for the computational details. We 
suppressed the effect of neutral-hydrogen collisions to a 
negligible level by multiplying the classical value based on 
Drawin's cross section (cf. Steenbock & Holweger 1984) 
by a factor of 10 -4 according to the conclusion of that 
paper. 

The derivation of the K abundance from a given EW 
was done as in Paper I with Kurucz's (1993) WIDTH9 
program, which had been considerably modified in vari- 
ous respects (especially to include the effect of departure 
from LTE). According to Takeda et al. (1996), the van 
der Waals damping parameter (Cg) was increased by ap- 
plying a correction of AlogCe = 4-1.00 to the conven- 
tional Unsold's (1955) formula value (corresponding to 
using T 6 = 2.5rU' lsold ). We adopted the log gf values of 
-0.13 (7665 line) and -0.17 (7699 line), and the radi- 



ation damping constant of r ra( j = 0.38 x 10 8 s 1 , which 
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Fig. 3. Spectra of 17 stars in the neighborhood of the K I 
7665 line (left panel) and the K I 7699 line (right panel). The 
seven spectra for the two standard stars and five M 13 stars in 
the left panel (drawn in blue lines) are the ones corrected for 
the telluric lines (cf. figure 2). Each spectrum (normalized 
with respect to the continuum) is vertically offset by 0.5 rel- 
ative to the adjacent one. An appropriate Doppler correction 
has been applied to each spectrum so that the positions of 
stellar spectral lines are located at the corresponding labora- 
tory wavelengths. 

were taken from Kurucz and Bell's (1995) compilation. 

Actually, we prepared two sets of non-LTE departure 
coefficients (grid of 90 models) corresponding to two dif- 
ferent choices of input K abundances ([K/Fc] = 0.0 and 
[K/Fe] = +0.5); considering that the resulting [K/Fe] is 
likely to be encompassed by these two values. And two 
kinds of non-LTE abundances (A^° and A^ 5 ) were ob- 
tained for a given EW for each of the two sets. Then, these 
and A 1 ^ 5 were interpolated (or extrapolated) so that 
the final non-LTE solution (Ak) and the used departure 
coefficient becomes consistent with each other (cf. sub- 
section 4.2 in Takeda & Takada-Hidai 1994). In addition, 
we also derived the LTE abundance A^ TE , from which the 
non-LTE correction was derived as A NLTE = Ak — A K TE . 
Finally, the differential K abundance relative to the Sun 
was computed as [K/H] = Ak — 5.12, where Anders and 
Grevesse's (1989) Ak,q value of 5.12 was adopted as the 
solar potassium abundance (which is also consistent with 
the result of Takeda ct al. 1996). The resulting values of 
[K/H] and A NLTE for each line along with the [K/Fe] ra- 
tio (= [K/H] - [Fe/H]; [K/H] is the average of two lines) 



are summarized in table 1. Regarding the typical errors 
in the K abundances caused by uncertainties in the atmo- 
spheric parameters or the damping constant, table 1 in 
Paper I (especially for the cases of metal-poor K giants 
such as HD 122563, HD 165195, and HD 221170) may be 
informative. 

5. Discussion 

5.1. Verifying Model Atmosphere Parameters 

It may be worth comparing our spectroscopically de- 
termined atmospheric parameters with those by another 
method occasionally used for globular cluster stars; i.e., 
the photometric and evolutionary logg ovo , where T e g 
is derived from colors and g is evaluated from L (the stellar 
luminosity estimated from the apparent magnitude and 
the distance), T e ff, and M (the stellar mass often assumed 
to be 0.8 Mq for globular cluster giants) by the relation 
of g/g Q = (M/M )(L/L )- 1 (T cff /r cff , ) 4 . 

The basic data of V magnitude and B — V color were 
taken from Ivans et al. (1999; M 4), Smith and Briley 
(2006; M 13), and Sneden et al. (1997; M 15); the color 
excess Eb-v and the distance d (assumed to be the same 
for all targets in a cluser) were taken from the on-line 
database 4 elaborated by the Padova group. For p Boo 
and a Boo, Eb-v was assumed to be zero, and the other 
data were taken from the SIMBAD database. We first 
derived T e ff from the dereddened (B — V)o color by using 
Alonso, Arribas, and Martinez-Roger's (1999) equation 
(4) in their table 2, where 0.0 (p Boo), -0.6 (a Boo), -1.2 
(M 4 stars), -1.7 (M 13 stars), and -2.5 (M 15 stars) were 
assumed for [Fc/H]. Then, since the stellar luminosity L 
can be obtained from the extinction-corrected Vq magni- 
tude (V — 3.1Eb~v), the distance d, and the bolometric 
correction evaluated by Alonso, Arribas, and Martinez- 
Roger's (1999) equations (17) and (18), we can derive 
logg cvo from the above-mentioned relation by assuming 
M = Q.8M Q . The resulting T e p ff ho and log# cvo (along with 
the corresponding [Fe/H] values when these parameters 
were used; cf. table 3) are summarized in table 2, where 
the published values taken from various literature (since 
1980's) are also presented for comparison. Besides, table 

3 gives the differences of these parameters relative to the 
adopted spectroscopic ones, and the resulting variations 
in the abundances of Fe as well as K caused by these pa- 
rameter changes. 

We can see the following tendency from an inspection 
of table 3: 

In the mildly metal-poor (M 4) or near-solar metallic- 
ity (a Boo, p Boo) region, we can not see any remark- 
ably systematic difference in §T e s(= — T^fT) or m 
<51og<?(= log 5 ovo - logg spc ); i.e., ST eS - ± 100 K and 
<51og 5 ~± 0.2-0.3 dex. 
- However, in the very metal-poor regime of [Fe/H] < 
-1.5 (M 13 and M 15), the inequality trend of ST eS > 
as well as <51og<? > manifestly appears, especially for the 
most metal-poor case of M 15 where X^ /log<? spe is appre- 

4 (http:/ /dipastro. pd.astro.it/globulars/). 
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Fig. 4. Spectra of the 17 program stars in the 6559—6567 A 
region comprising Ha. Each spectrum is normalized with re- 
spect to the continuum (horizontally-drawn red dashed line) 
and is vertically shifted relative to the adjacent one by an ap- 
propriate offset (1.0 or 0.5). The wavelength scale is adjusted 
to the same system as in figure 3 where the stellar metallic 
lines locate at the laboratory wavelengths. The position of 
the Ha line center (6562.8 A) in this reference frame is indi- 
cated by the vertically-drawn blue dashed line. The redward 
emission for III- 73 (M 13) is so strong (peaked at ~ 1.6) that 
it is crossing the profile of 111-59 just above. 

ciably lower than T c P g°/log 5 cvo by ~ 100-200 K/~ 1 dex 
(even \ogg spc < cases are seen, which is hardly possible 
for hydrostatic stellar atmospheres). 

Given the existence of such an appreciable discrepancy 
we tend to wonder whether our use of spectroscopic pa- 
rameters is really justified. Which sets (T^ c /\ogg spc vs. 

^eff °/l°g.9 cvo ) should be preferably used for abundance 
determinations? Before going into this, we point out that 
that our targets of intrinsically bright red giants are on 
or near to the AGB tip where stars (in their late stage of 
evolution) are likely to have extended envelopes and show 
significant mass loss as well as time variability. Actually, 
as demonstrated in figure 4, not a few stars show promi- 
nent emissions and blue-shifted cores in Ha (especially, 
emissions are seen in all five M 15 stars), indicating the 
existence of such active phenomena. 

In such cases, we would state that the spectroscopic 
method is more advantageous: 

(1) If stars are unstable and variable in time, parame- 
ters based on observational data (colors, magnitudes, etc.) 
simply taken from catalogues are not reliable any more. 
Meanwhile, spectroscopic parameters are established from 
the spectrum itself, from which abundances are derived. 

(2) We must recall that most stellar model atmospheres 
widely used are based on the assumption of hydrostatic 
equilibrium and the plane-parallel approximation. In case 
where these assumptions failed, it would not be sensible 
any more to assign parameters derived from the true fun- 
damental stellar quantities (e.g., logg 6 ™ derived from re- 
alistic L, M, ...). On the other hand, it is still possible to 
apply conventional model atmospheres, if we could choose 
their parameters so carefully as to reproduce the physical 
condition of the real atmosphere (while regarding that T e g 
and logg are not so much real physical quantities as rather 
adjustable parameters). Our T e s | c and log<7 spc should thus 
be interpreted in this sense. 



Consequently, according to our opinion, the appreciably 
low T*g C and low log<7 spo in M 13 and M 15 are nothing 
but a manifestation of the existence of extended cooler 
region of lower density affecting the formation of spectral 
lines. As far as we invoke classical plane-parallel model 
atmospheres, such low parameter values must have been 
assigned to reproduce the relevant atmospheric condition. 
Since this kind of discordance is pronounced in M 13 

([Fe/H] 1.7) and especially in M 15 ([Fe/H] 2.5) 

while not in M 4 ([Fe/H] ~ — 1.1), this effect is considered 
to be metallicity-dependent, which may presumably be re- 
lated to the observational fact pointed out by Meszaros, 
Dupree, and Szentgyorgyi (2008) that the outflow veloci- 
ties of red giants are higher in metal-poor clusters (M 15) 
than in metal-rich clusters (M 4). To sum up, we consider 
that the use of and low logg spc is surely reasonable 
and preferable to other choices. 5 

5.2. Homogeneity of Fe and K Abundances 

The results of [Fe/H], [K/H] 7665 , and [K/H] 769 9 (cf. ta- 
ble 1) for each star are graphically depicted in figure 5, 
from which we can read the following characteristics: 

— The [Fe/H] values are considered to be nearly uniform 
in all three clusters (M 4, M 13, and M 15) to within a 
precision of ~ 0.05-0.1 dex. Actually, the mean ([Fe/H]) 
(and the standard deviation ±er) averaged over 5 stars are 
-1.14(±0.04), -1.65(±0.09), and -2.48(±0.07) for M 4, 
M 13, and M 15, respectively. 6 Therefore, we can state 
that the any of our M 4/M 13/M 15 targets have essen- 
tially the same metallicity within the cluster. 

- Regarding [K/H], for which values derived from K I 
7655 and 7699 lines are consistent with each other, we 
notice that two stars (M 4/L2406 and M 13/111-73) are 
evidently anomalous because they show considerably dis- 

5 As an alternative interpretation, it may be possible to consider 
that the large discrepancy in logg (systematically low logg spe ) 
seen very metal-poor red giants is due to the non-LTE overion- 
ization effect (Fe I lines are weakened compared to the case of 
LTE, while Fe II lines are practically unaffected), which would 
naturally yield an underestimation of log</ spe if LTE is assumed. 
We can not exclude this possibility, as Ruland et al. (1980) 
once reported a sign of non-LTE overionization in low-excitation 
Fe I lines (while high-excitation lines are comparatively inert) 
in their analysis of the early-K giant Pollux (T e £f ~ 4800 K, 
logg ~ 2.2, and the near-solar metallicity). On the other hand, 
theoretical investigations on the non-LTE effect in the forma- 
tion of Fe lines in red giants are difficult, because the results 
are sensitively influenced by uncertainties in computational de- 
tails (e.g., the treatment of UV photoionizing radiation field, 
collisional cross section with neutral hydrogen atoms; see, e.g., 
Steenbock 1985 or Takeda 1991). As far as the present case 
is concerned, however, we consider that the possibility of an 
appreciable non-LTE overionization in Fe I such as claimed by 
Ruland et al. (1980) is rather unlikely. That is, since the charac- 
teristic feature they found was "low-excitation Fe I lines tend to 
yield lower LTE abundances than high-excitation Fe I lines" , the 
LTE excitation temperature would appear higher than usual, by 
which an overestimation of T^ c should eventually result. This 
is, however, just the opposite to what we found in M 13 or M 15 
stars. 

6 For reference, a simple averaging of the literature data given in 
table 2 (excluding our results) gives the mean cluster metallicties 
of 1.2 (M 4), ~ -1.5 (M 13), and 2.3 (M 15). 
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Fig. 5. Graphical plots of the Fe and K abundances (relative 
to the Sun) derived for each of the 17 stars (also given in 
table 1). The crosses (connected by dashed lines for those 
5 stars belonging to the same cluster), filled triangles, and 
open squares indicate [Fe/H], [K/H] 7 6g5 (from K I 7665 line), 
and [K/H] 76 99 (from K I 7699), respectively. Plots for M 4, 
M 13, and M 15 stars are represented in red, green, and blue, 
respectively. 

crcpant (i.e., larger) [K/H] compared to other members. 
However, when these two stars are excluded, we can con- 
firm a reasonable uniformity of [K/H] within a cluster; i.e., 
<[K/H] 7665 )/([K/H] 7699 ) (±a 7665 /±a 7699 ) is -1.01/- 1.03 
(±0.04/±0.04) for M 4 (4 stars), -1.51/— (±0.20/—) for 
M 13 (4 stars), and -2.26/ -2.23 (±0. 12/ ±0.09) for M 15 
(5 stars). The larger dispersion (0.20) for the case of M 13 
compared to other two clusters (~ 0.05-0.1) may be un- 
derstood as due to the availability of only one K I 7665 
line which is severely contaminated by telluric lines (cf. 
figure 2). 

How should the discrepant nature of these two outliers 
be interpreted? Actually, we have a good reason to be- 
lieve that this is nothing but a superficial effect caused 
by an inadequate treatment in the abundance determi- 
nation (i.e., not the real abundance anomaly). We point 
out that the spectra of these M 4/L2406 and M 13/111-73 
have markedly broad line widths (figure 3) and promi- 
nently strong Ha emission components (figure 4) than 
the remaining stars, which suggests a significant increase 
in the activity or turbulent velocity fields in the upper 
atmosphere (where the core of K I lines form) presum- 
ably related to dynamical phenomena often seen in AGB 
stars. We thus suspect the existence of considerably 
dcpth-dcpcndcnt turbulent velocity dispersion which in- 



creases with height, 7 which may have contributed an ad- 
ditional increase to the strength of high- forming saturated 
K I lines. In such a case, an application of the v t value 
determined from Fc lines with EW < 100 mA would be 
no more adequate for the K I lines with EW of a few 
x 100 mA (especially intensified by the enhanced tur- 
bulence in the upper atmosphere), for which a somewhat 
larger vt value should have been more relevant. Hence, we 
believe that the reason for the anomalously large [K/H] 
obtained for these two stars is the assignment of inappro- 
priately small f t values, which must have resulted in an 
overestimation of K abundances from such strongly satu- 
rated K i 7665/7699 lines. 

Consequently, regarding that M 4/L2406 and M 13/111- 
73 correspond to exceptional cases of peculiar velocity 
fields in the upper atmospheric layer, we conclude that 
the abundance uniformity essentially holds for [K/H] (as 
well as for [Fe/H]) within M 4, M 13, and M 15. 

5.3. Behavior of [K/Fe] over Wide Metallicities 

Now that the homogeneity of K (and Fe) abundance 
within each globular cluster has been confirmed, we can 
regard based on the argument in section 1 that the ob- 
served [K/Fe] ratios (which naturally turn out also uni- 
form) of the cluster stars retain the original composition 
of the halo gas, from which they formed. 

According to table 1 (where the mean of [K/Fe] 7 665 
and [K/Fe] 7 6 99 is given), the averaged cluster val- 
ues of ([K/Fe]) (and the standard deviation ±er) are 
+0.10(±0.04) (M 4, [Fe/H] ~ -1.1, 4 stars, L2406 ex- 
cluded), ±0.12(±0.16) (M 13, [Fe/H] ~ -1.7, 4 stars, III- 
73 excluded), and +0.24(±0.05) (M 15, [Fe/H] ~ -2.5, 
5 stars). We can state from these results that [K/Fe] is 
marginally supersolar at ~ +0.1-0.3 and [K/Fc] tends to 
gradually increase with a decrease in the mctallicity, which 
means a rather tight and monotonic [K/Fe] vs. [Fe/H] re- 
lation (with almost no significant diversity) in the range 
of -2.5 < [Fc/H] < -1. These [K/Fe] vs. [Fe/H] rela- 
tions for 15 stars of M 4, M 13, and M 15 (along with 
p Boo and a Boo) are plotted in figure 6, where the re- 
sults for nearby F-G disk dwarfs (— 1 < [Fc/H] < 0; taken 
from Paper I) and those for extremely metal-poor stars 
(—4 < [Fe/H] < —2.5; obtained by reanalyzing the EW 
data of Cayrel et al. 2004; cf. Appendix) are also shown 
for comparison. 

As can be clearly seen from this figure, our [K/Fe] 
results in this study for the cluster stars (excluding 
M 4/L2406 and M 13/111-73) and the reference objects 
(p Boo and a Boo) over — 2.5 ^ [Fc/H] < connect with 
those of extremely mctal-dcficicnt stars as well as disk 
stars quite satisfactorily. We may thus conclude that 
[K/Fe] behaves itself in a rather clear manner over a 
wide [Fe/H] range of ~ 4 dex from near-solar to ultra- 
low metallicity. That is, the gradual increase of [K/Fe] 

7 More generally speaking, such an increasing tendency of the 
turbulent velocity field with height is a phenomenon occasionally 
seen in low-gravity stars (though its degree is different from case 
to case); see, e.g., Takeda (1992) for a Boo (K giant) or Takeda 
and Takada-Hidai (1994; appendix B) for A— F supergiants. 
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Fig. 6. (a) The [K/Fe] vs. [Fe/H] plots constructed from 
the metallicity and the NLTE potassium abundances given in 
table 1. The results for the 17 program stars are shown in 
large open symbols (diamonds for p Boo and a Boo, squares 
for M 4 stars, circles for M 13 stars, and triangles for M 15 
stars). The two stars (M 13 III- 73 and M 4 L2406) show- 
ing exceptionally large deviations from the main trend are 
also indicated. Meanwhile, the smaller filled symbols repre- 
sent the results for disk stars (circles, taken from Paper I) 
and extremely metal poor stars (inverse triangles, reanalysis 
of Cayrel et al.'s 2004 data, see the Appendix) presented for 
comparison. Lines show the three kinds of theoretical pre- 
dictions (see section 1 of Paper I for more details): Dashed 
line (marked "S") — Samland (1998; taken from his figure 
14), dash-dotted line (marked "T") — Timmes, Woosley, and 
Weaver (1995; the top curve in their figure 24), solid line 
(marked "G") — Goswami and Prantzos (2000; "Case B" in 
their figure 7). (b) The NLTE corrections for K abundance 
determinations (A NLTB is the average of Ai^Tp and 
in case that both lines are available), plotted against [Fe/H]. 
The same meanings of the symbols as in (a). 

with a decrease in [Fc/H] (d[K/Fe]/<f[Fe/H] 0.3) seen 

in disk stars of — 1 < [Fe/H] < (Paper I) continues fur- 
ther down at [Fe/H] < — 1 with a slightly decreased slope 
until [Fc/H] ~ —3, where [K/Fe] appears to attain a broad 
maximum of ~ +0.3 (followed by a sign of weak downturn 
of [K/Fc] at [Fc/H] < -3). 

Accordingly, we do not regard the considerably large 
diversity of [K/Fe] previously reported for halo stars of 
[Fe/H] - -2 (Zhang et al. 2005, Gratton & Snedcn 1987a) 
to be real, which we suspect to be due to some improper 
treatment in abundance determinations (e.g., EW errors 
caused by unsuccessful removal of telluric lines, or use of 
an inadequate microturbulence such as we encountered in 
the present cases of M 4/L2406 and M 13/111-73). 

This observational fact of mildly supersolar [K/Fe] over 
a wide metallicity range (with [K/Fe] ~ at near-solar 
metallicity and a broad peak at [Fe/H] ~ —3) has to 
be explained by theoretical calculations. However, the- 
oreticians did not find it easy to reproduce the trend of 
[K/Fe] > 0, especially at the very metal-poor regime, be- 
cause the K yield is metallicity-dcpcndent (as a charac- 
teristic of odd-Z element) and becomes too small there 
to bring [K/Fe] above zero. Therefore, comparatively 
plain calculations predicted only subsolar (or near-solar 



at most) [K/Fe] ratios for metal-poor stars. For exam- 
ple, Timmes, Woosley, and Weaver's (1995) calculations 
resulted in a decrease of [K/Fe] 0) with a lowering of 
[Fe/H] as shown figure 24 of their paper (also plotted in 
figure 6 as curve "T" ) , which markedly contradicts the 
observed tendency. Therefore, if the supersolar [K/Fe] is 
to be reenacted even at the metal-poor regime, one had 
to contrive ways to raise it by some means, such as empir- 
ically adjusting the K yield (Samland 1998) or using an 
IMF of larger exponent (cf. Goswami & Prantzos 2000; 
putting larger weight to comparatively less massive SNcII 
of lower gravitational potential), as described in section 1 
of Paper I, by which a trend more or less similar to the 
observed one could somehow result (cf. figure 6 plotted as 
curves "S" and "G", respectively). At any rate, it is evi- 
dent that our understanding on the synthesis mechanism 
of K and its evolution in the Galaxy is still considerably in- 
sufficient (see also the Appendix for the K yield problem in 
the extremely metal-poor regime). Now that we have be- 
come more confident than before about the observational 
[K/Fc] vs. [Fe/H] trend over a wide metallicity range, we 
would encourage theoreticians to revisit this problem, so 
that a better match between theory and observation may 
be established in a reasonable manner. 

6. Conclusion 

Given the confusing situation regarding the [K/Fc] ratio 
of metal-poor stars in the halo, where some work suggests 
a fairly tight tendency (such like disk stars) while others 
report a considerably large diversity amounting to ~ 1 dex, 
we carried out a spectroscopic study on 15 red giants of 
three mildly to very metal-poor globular clusters (M 4, 
M 13, and M 15) along with two reference stars (p Boo and 
a Boo) based on the high-dispersion spectra obtained with 
Subaru/HDS, with a purpose of clarifying the behavior of 
[K/Fe] at the relevant metallicity range of —2.5 < [Fe/H] < 
-1. 

The atmospheric parameters (T e g, \ogg, Vt, and [Fe/H]) 
were spectroscopically determined by using Fe I and Fe II 
lines, and the abundance of K was derived from the K I 
resonance lines at 7664.91 and 7698.97 A while taking into 
the non-LTE effect. 

We confirmed that [K/H] (as well as [Fe/H]) is almost 
homogeneous within each of the three clusters to a pre- 
cision of ^ 0.1 dex, though superficially large deviations 
are exceptionally seen for two peculiar stars which show 
signs of considerably increased turbulence in the upper 
atmosphere. 

The [K/Fe] ratios of these cluster stars turned out 
mildly supersolar by a few tenths of dex, tending to grad- 
ually increase from ~ +0.1-0.2 at [Fe/H] ~ -1 to ~ +0.3 
at [Fe/H] ~ —2.5, which is a fairly tight and clean trend. 
We thus consider that the previously reported large diver- 
sity of [K/Fe] in halo stars is not real, which we suspect 
to be due to some improper treatment in the analysis. 

This result is quite consistent (i.e., smoothly connect- 
ing) with the [K/Fe] trend of disk stars (—1 < [Fc/H]) and 
that of extremely metal-poor stars (— 4< [Fc/H] < —2.5). 
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That is, [K/Fe] appears to continue a gradual increase 
from [Fc/H] ~ toward a lower metallicity regime down 
to [Fe/H] ~ —3, where a broad maximum of [K/Fe] ~ 
+0.3-0.4 is attained, possibly followed by a slight down- 
turn toward a further lower metallicity at [Fe/H] ^ — 3. 

This investigation is based on the data obtained by the 
observation with the Subaru Telescope, which was car- 
ried out during the practical training of observational as- 
tronomy for graduate students as a courscwork of The 
Graduate University for Advanced Studies (SOKENDAI). 
We express our hcartful thanks to S. S. Hayashi, R. 
Furuya, and A. Tajitsu for their continuous support and 
encouragement, as well as to C.-H. Peng and S. Honda for 
their collaboration in the observation. 

Helpful comments by N. Tominaga and N. Iwamoto 
from the theoretical side on the first version of the 
manuscript are also acknowledged. 

This research made use of the SIMBAD database oper- 
ated by CDS, Strasbourg, France. 

Appendix. Reanalysis of Cayrel et al.'s (2004) 
Data 

In the extensive spectroscopic study on the sample of 
35 extremely metal-poor stars toward clarifying the abun- 
dance patterns of 17 elements from C to Zn recently 
carried out by Cayrel et al. (2004), the abundances of 
K were also determined by using K I 7665/7699 lines. 
However, they applied a constant non-LTE correction of 
A NLTE = —0.35 dex to all stars, which they adopted as a 
rough estimate based on Ivanova and Shimanskii's (2000) 
calculations. Because of this apparently imperfect treat- 
ment of the non-LTE effect, one can not be sure whether 
the behavior of [K/Fe] they obtained (i.e., near-solar or 
slightly supersolar with a decreasing tendency toward a 

lower metallicity from [K/Fe] ~ +0.2 at [Fe/H] 2.5 to 

[K/Fe] ~ at [Fe/H] 4; cf. their figure 9) is reliable 

or not. 

We decided, therefore, to reanalyze their EW data of 
K I 7665/7699 lines with the atmospheric parameters they 
used, while properly taking into account the non-LTE ef- 
fect. Regarding the non-LTE departure coefficients to be 
included in line-formation calculations for abundance de- 
termination, we used the already calculated results for a 
grid of models (T off from 4500 K to 6500 K, log 5 from 
1.0 to 5.0, v t from 1 to 3 km s _1 , and [Fe/H] from 0.0 to 
—3.0), which were originally computed for constructing 
extensive tables of non-LTE corrections 8 as explained in 
the Appendix of Paper I. Besides, we newly performed 
calculations corresponding to the case of [Fe/H] = —4 
to be combined with the previous grids, so that we can 
handle any of the relevant 31 stars (for which Cayrel et 

8 The anonymous ftp site to access these electronic 
tables described in the Appendix of Paper I is 
now not available any more. Instead, the same 

data materials are placed at the following web site: 
(http: / /optik2. mtk.nao.ac.jp/~takeda/potassium_nonltc/). 



al. measured -EW7665 and/or EWjqg^) in the parame- 
ter range of 4500 K < T cff < 5300 K, 0.7 < logs < 2.7, 
1.2kms- 1 <-y t <2.2kms- 1 , and -4.0 < [Fc/H] <-2.1 by 
interpolation. Following the same manner as described in 
subsection 4.2, we determined [K/H], A NLTE , and [K/Fe] 
for each star, as given in table 4. The resulting values of 
[K/Fe] (as well as A NLTE ) are plotted against [Fe/H] in 
figure 6 (filled inverse triangles). 

As we can see from figure 6b, the extents of the non- 
LTE corrections ( | A NLTE | ) systematically decrease with a 
lowering of the metallicity; such as from A NLTE ~ —0.3 (at 
[Fe/H] ~ -2.5) to A NLTE ~ -0.2 (at [Fe/H] ~ -4), which 
means that Cayrel et al.'s (2004) use of A NLTE = -0.35 
must have somewhat overcorrected (i.e., underestimated) 
[K/Fc] by ~ 0.1-0.2 dex, such that the amount of overcor- 
rection progressively increasing with a decrease in [Fc/H]. 
Accordingly, regarding the [K/Fe] vs. [Fe/H] relation of 
extremely metal-poor stars, we consider that Cayrel et 
al.'s (2004) result ([K/Fe] ratio is near-solar or slightly 
above solar (~ +0.1-0.2) with an appreciable gradient 
of d[K/Fe]/d[Fe/H]) should be revised as concluded in 
subsection 5.3 (cf. figure 6a): [K/Fe] in this very metal- 
deficient regime still remains well supersolar at ~ +0.2- 
0.3, with a marginal sign of decline (downturn) toward a 
further lower [Fe/H]. 

The conclusion that [K/Fe] remains supersolar even at 
[Fe/H] ~ —4 may provide theoreticians with an important 
constraint on the synthesis mechanism of K from the ob- 
servational side, because considerable difficulties exist in 
theoretically reproducing the [K/Fe] ratio (cf. subsection 
5.3) also in this distinctly low metallicity regime, where 
the observed [K/Fe] may be considered to simply reflect 
the composition of the first super/hyper-nova ejecta. For 
example, Tominaga, Umeda, and Nomoto's (2007) nucle- 
osynthesis calculation based on population III supernova 
models predicted [K/Fe] ~ -1(±0.5) at -4< [Fc/H] < -3 
(cf. their figure 6), which apparently disagrees with the 
observational fact mentioned above. Hence, in order to 
resolve this situation, a new reassessment of assumptions 
or fundamental physics might as well be required. In this 
connection, Iwamoto et al.'s (2006) finding, that odd-Z 
elements such as K and Sc can be significantly overpro- 
duced if the proton-rich environment (large initial electron 
fraction Y e ) is realized, is quite interesting. Further inves- 
tigation following this line may be worth a try. 
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Table 1. Stellar parameters and potassium abundance results. 



Star 


T e g 


loggr 


Ut 


[Fc/H] 


£W 7 665 [K/H] 7665 


A N LIE 

7665 


-EW7699 


K/H] 7 699 


A NLTE 

7699 


[K/Fc] 




(K) (cm 


s" 1 )(k 


m 


') (dex) 


(mA) 


(dex) 


(dex) 


(mA) 


(dex) 


(dex) 


(dex) 


p Boo 


4363 


2.09 


1.29 


+0.05 


*287.7 


-0.19 


-0.08 


238.6 


-0.25 


-0.12 


-0.27 


a Boo 


4281 


1.72 


1.49 


-0.55 


284.5 


-0.45 


-0.12 


235.9 


-0.55 


-0.19 


+0.05 


M4 L2406 


4048 


0.61 


2.16 


-1.22 


364.3 


-0.39 


-0.09 


319.5 


-0.42 


-0.12 


+0.81 


M4 L2617 


4256 


1.48 


1.38 


-1.12 


223.1 


-0.99 


-0.21 


197.9 


-0.98 


-0.26 


+0.13 


M4 L3209 


4025 


1.13 


1.49 


-1.12 


255.4 


-0.98 


-0.13 


222.5 


-1.01 


-0.18 


+0.12 


M4 L3624 


4269 


1.47 


1.42 


-1.12 


216.8 


-1.07 


-0.24 


189.5 


-1.10 


-0.30 


+0.04 


M4 L4511 


4173 


1.20 


1.44 


-1.13 


230.8 


-0.99 


-0.18 


202.9 


-1.02 


-0.23 


+0.12 


M13 1-13 


4155 


0.73 


1.75 


-1.58 


*215.4 


-1.36 


-0.26 








+0.22 


M13 11-76 


4159 


0.27 


1.87 


-1.77 


*186.8 


-1.70 


-0.25 








+0.07 


M13 111-52 


4271 


0.80 


1.55 


-1.58 


"168.1 


-1.67 


-0.36 








-0.09 


M13 111-59 


4206 


0.54 


1.42 


-1.58 


* 190.7 


-1.31 


-0.26 








+0.27 


M13 111-73 


4145 


0.38 


1.93 


-1.73 


*252.8 


-1.03 


-0.23 








+0.70 


M15 K144 


4066 - 


-0.45 


1.46 


-2.51 


108.2 


-2.28 


-0.17 


90.9 


-2.24 


-0.18 


+0.25 


M15 K341 


4085 - 


-0.15 


2.10 


-2.51 


131.2 


-2.26 


-0.20 


102.0 


-2.28 


-0.19 


+0.24 


M15 K431 


4152 - 


-0.01 


1.49 


-2.47 


117.8 


-2.12 


-0.28 


89.4 


-2.21 


-0.25 


+0.31 


M15 K634 


4161 


0.23 


1.78 


-2.36 


126.4 


-2.19 


-0.27 


111.8 


-2.09 


-0.27 


+0.22 


M15 K825 


4012 - 


-0.22 


1.78 


-2.55 


112.9 


-2.44 


-0.13 


101.3 


-2.31 


-0.15 


+0.17 



Note. 

Columns 2-5 give the atmospheric parameters (the effective temperature, the surface gravity, the microturbulence, and the 
Fe abundance relative to the Sun) which were spectroscopically determined by using Fe I and Fe II lines and adopted in this 
study. The line equivalent width, the non-LTE potassium abundance relative to the Sun, and the relevant non-LTE correction 
are presented in columns 6-8 (K I 7665 line) and 9-11 (K I 7699 line). In the last column 12 is given the K-to-Fe logarithmic 
abundance ratio, [K/Fe] (= [K/H] — [Fc/H]), where [K/H] is the average of [K/H]7665 and [K/H]7699 in case both lines are 
available. 

* Corrected values for the imperfect removal of telluric lines; AEW ICB id corrections of 38.8 mA (p Boo), 30.9 mA (1-13), 37.9 mA 
(11-76), 47.9 mA (111-52), 72.9 mA (111-59), and 22.8 mA (III- 73) have been subtracted from the directly measured EW (cf. 
footnote 3 in subsection 4.1). 
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Table 2. Comparison of the atmospheric parameters and metallicity with the literature values. 



Star 'J^ff 




Vt 


[Fe/Hj 


Kclcrcnccs 


(K) 


(cm s~ 1 ) 


(km s" 1 ) 


(dcx) 




p Boo 4363 


2.09 


1.29 


+0.05 


This study (adopted, spectroscopic parameters) 


1 2 > 2 






0.00 


This study (T^n° and log(7 GV °^ f°r reference) 


4260 


2.22 


2.10 


— 0.17 


MC W90 


a Boo 4281 


1.72 


1.49 


-0.55 


this study (adopted, spectroscopic parameters) 


4303 


1.56 




-0.65 


This study (T^° and log(/ Gvo1 for reference) 


4490 


2.01 


1.80 


-0.56 


LAM81 


112", 


1.06 


2.50 


-0.48 


GRA82 


4350 


1.60 


1.70 


-0.58 


BEL85 


4330 


1.50 


1.50 


-0.38 


GRA86 


4400 


1.70 


2.30 


-0.55 


KYR86 


4250 


1.70 


2.40 


-0.60 


LEE87 


4375 


1.97 


1.80 


-0.42 


EDV88 


4300 


2.00 


1.50 


-0.69 


FER90 


4280 


2.19 


2.30 


-0.60 


MCW90 


4330 


2.10 


1.60 


-0.58 


BR092 


4280 


1.30 


1.40 


-0.54 


MCW94 


4300 


1.50 


1.70 


-0.47 


SNE94 


4300 


1.50 


1.70 


-0.51 


HIL97 


4250 


1.30 


1.70 


-0.68 


GON98 


4345 


2.05 


1.50 


-0.37 


THE99 


4300 


1.50 


1.70 


-0.63 


TOM99 


4300 


1.50 


1.72 


-0.49 


CAR00 



M4 L2-1IH) 



4048 
4125 
4100 
4150 



0.61 
0.76 
0.45 
0.15 



2.16 



2.45 
2.20 



-1.23 This study (T^?° and logg Gvo1 for reference) 



-1.20 
-1.30 



IVA99 
YON08 



M4 L2617 



4256 
4238 
4200 
4275 



1.48 
1.24 
0.95 
1.25 



1.38 



1.55 
1.65 



-1.12 This study (adopted, spectroscopic parameters) 

— 1.21 This study (T p ^° and log(/ Gvo1 for reference) 

-1.17 IVA99 

-1.20 YON08 



M4 L3209 4025 
4162 
3975 

4075 

M4 L3624 4269 
4290 
4225 
1225 



1.13 
0.83 
0.60 
0.75 

irw 

1.26 
1.10 
1.05 



1.49 —1.12 This study (adopted, spectroscopic parameters) 



and log(/ GVO for reference) 



-1.33 This study (T P j° 
1.75 -1.20 IVA99 

1.95 -1.25 YON08 

This study (adopted, spectroscopic parameters) 



1.45 
1.60 



-1.12 

-1.22 This study (T P j° 

-1.16 IVA99 

-1.29 YON08 



and lo 



for reference) 



M4 L4511 



4173 
4251 
4150 
4150 



1.20 
1.16 
1.10 
1.05 



1.44 



1.55 
1.70 



-1.13 
-1.21 
-1.19 
-1.22 



This study (adopted, spectroscopic parameters) 



This study (T 
IVA99 
YON08 



pho 



and logg GVO for reference) 



M13 1-13 4155 
4200 

4290 

M13 11-76 4159 
4212 
4350 
4300 



0.73 
0.91 
1.00 
0.27 
0.91 
1.00 
0.85 



1.75 —1.58 This study (adopted, spectroscopic parameters) 

■■■ —1.57 This study (T p f ^° and logg Gvo1 for reference) 

2.00 -1.46 * 

"T8? -1.77 



KRA93, KRA92, SHE96 

This study (adopted, spectroscopic parameters) 



— 1.65 This study (T^° and logg Gvo1 for reference) 
2.00 -1.49 KRA93. KRA92. SHE96. SMI96 
1.95 -1.53 COH05 



M13 111-52 4271 0.80 1.55 -1.58 This study (adopted, spectroscopic parameters) 

4248 1.00 ■■■ -1.56 This study (T P j° and log g° vo1 for reference) 

4335 1.00 2.00 -1.50 KRA93, KRA92, SHE96 



M13 111-59 



M13 111-73 



4206 
4260 
4360 
4145 
4177 
4300 



0.54 
0.97 
1.10 



0.38 
0.81 
0.85 



1.42 —1.58 This study (adopted, spectroscopic parameters) 

--- —1.50 This study (T p f ° and log(/ Gvo1 for reference) 

1.75 -1.45 KRA93, KRA92, SHE96, SMI96 

1.93 —1.73 This study (adopted, spectroscopic parameters) 

■■■ —1.66 This study (T p |? and log(/ Gvo1 for reference) 

2.25 -1.51 KRA93, KRA92, SHE96, SMI96 



M15 K144 



4066 
4242 
4390 
4460 
4425 



-0.45 
0.74 
0.90 
0.95 
0.75 



1.46 

2.00 
2.00 
2.10 



-2.51 
-2.25 
-2.31 
-2.21 
-2.34 



This study (adopted, spectroscopic parameters) 



ph., 



This study (T p „ 
SNE91 
SNE91 
SNE97 



and log(/ GVO for reference) 



M15 K341 



4085 
4210 
4275 



-0.15 
0.61 
0.45 



2.10 



2.00 



-2.51 
-2.35 
-2.34 



This study (adopted, spectroscopic parameters) 



pho 



This study (T p " 
SNE97 



:id logp GVO for reference) 



M15 K431 4152 


-0.01 


1.49 


-2.47 


This study (adopted, spectroscopic parameters) 


4284 


0.75 




-2.30 


This study (T p ° and logg Gvo1 for reference) 


4430 


0.90 


2.00 


-2.28 


SNE91 


4500 


1.00 


2.00 


-2.18 


SNE91 


4375 


0.50 


2.30 


-2.43 


SNE97 



M15 K634 



4161 
4242 
4225 
4225 



0.23 
0.64 
0.30 
0.60 



1.78 



1.85 
2.05 



npho 



and logg Gvo1 for reference) 



This study 
SNE97 

OTS06 

This study (adopted, spectroscopic parameters) 



and logp GVO for reference) 



-2.36 This study (adopted, spectroscopic parameters) 
-2.27 
-2.34 
-2.30 

M15 K825 4012 -0.22 1778 -2.55 

4242 0.63 ■ ■ ■ -2.30 This study (T P j° 

4275 0.65 1.75 -2.42 SNE97 ° 

Note. In case that two kinds of logg values (g spc and g cvo ) arc given for the same star (e.g., IVA99 for M 4 stars 
or SNE97 for M 15 stars), we prefereably adopted the spectroscopic logg spe . 

Key to the references: BEL85 — Bell, Edvardsson, and Gustafsson (1985); BR092 — Brown & Wallerstein (1992); 
CAR00 — Carr, Sellgren, and Balachandran (2000); COH05 — Cohen & Melendez (2005); EDV88 — Edvardsson 
(1988); FER90 — Fernandez- Villacanas, Rego, and Cornide (1990); GON98 — Gonzalez & Wallerstein (1998); 
GRA82 — Gratton et al. (1982); GRA86 — Gratton & Ortolani (1986); HIL97 — Hill (1997); IVA99 — Ivans et al. 
(1999); KRA92 — Kraft et al. (1992); KRA93 — Kraft et al. (1993); KYR86 — Kyrolainen et al. (1986); LAM81 

- Lambert & Ries (1981); LEE87 — Leep, Wallerstein, and Oke (1987); MCW90 — McWilliam (1990); MCW94 

- McWilliam & Rich (1994); OTS06 — Otsuki et al. (2006); SHE96 — Shetrone (1996); SMI96 — Smith et al. 
(1996); SNE91 — Sneden et al. (1991); SNE94 — Sneden et al. (1994); SNE97 — Sneden et al. (1997); THE99 - 
Thevenin & Tdiart 119991: TOM99 — Tomkin Lamhert 11999): YON08 — Yonsr et al. 120081. 
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Table 3. Abundance variations in case that the photomeric T e g and the evolutionary logg are used. 



Star 




dlogg 


d[Fc i/HJ d[Fc n/HJ 


d[Fe/H]av 


d[K/Hj 7 6C5 


d[K/Hj 7 699 


d[K/HJ av 


i[K/Fe] av 




(K) 


(cm s _1 ) 


(dex) 


(dex) 


(dex) 


(dex) 


(dex) 


(dex) 


(dex) 


p Boo 


-111 


-0.34 


-0.06 


-0.04 


-0.05 


-0.05 


-0.06 


-0.06 


0.00 


a Boo 


+22 


-0.16 


-0.04 


-0.15 


-0.10 


+0.07 


+0.05 


+0.06 


+0.16 


M4 L2406 


+77 


+0.15 


+0.03 


—0.04 


0.00 


+0.06 


+0.07 


+0.06 


+0.07 


M4 L2617 


-18 


-0.24 


—0.05 


—0.13 


-0.09 


+0.03 


+0.02 


+0.03 


+0.12 


M4 L3209 


+137 


-0.30 


-0^02 


-0^40 


-0.21 


+0.25 


+0.22 


+0.24 


+0.44 


M4 L3624 


+21 


-0.21 


-0.02 


-0.17 


-0.10 


+0.07 


+0.05 


+0.06 


+0.16 


M4 L4511 


+78 


-0.04 


+0.01 


-0.17 


-0.08 


+0.11 


+0.11 


+0.11 


+0.19 


M13 1-13 


+45 


+0.18 


+0.02 


-0.02 


0.00 


+0.05 




+0.05 


+0.05 


M13 11-76 


+53 


+0.64 


+0.01 


+0.23 


+0.12 


-0.02 




-0.02 


-0.14 


M13 111-52 


-23 


+0.20 


-0.03 


+0.07 


+0.02 


-0.03 




-0.03 


-0.05 


M13 111-59 


+54 


+0.43 


+0.05 


+0.11 


+0.08 


+0.03 




+0.03 


-0.05 


M13 III— 73 


+32 


+0.43 


0.00 


+0.15 


+0.07 


-0.04 




-0.04 


-0.11 


M15 K144 


+ 176 


+1.19 


+0.16 


+0.37 


+0.26 


+0.08 


+0.10 


+0.09 


-0.17 


M15 K341 


+ 125 


+0.76 


+0.11 


+0.23 


+0.17 


+0.06 


+0.07 


+0.07 


-0.10 


M15 K431 


+ 132 


+0.76 


+0.13 


+0.22 


+0.17 


+0.08 


+0.09 


+0.09 


-0.09 


M15 K634 


+81 


+0.41 


+0.10 


+0.09 


+0.10 


+0.10 


+0.11 


+0.11 


+0.01 


M15 K825 


+230 


+0.85 


+0.28 


+0.24 


+0.26 


+0.20 


+0.21 


+0.21 


-0.05 



Note. 

Given are the abundance variations by using the photometric T e g and evolutionary logp instead of the spectroscopically 
determined "standard" parameters adopted in this study (cf. table 1). (Note that the same EW data for Fe as well as K lines 
and the same microturbulence were used in this test calculation.) Columns 2 and 3 give the relevant changes in T e g (T^° — Tj2 e ) 
and logg (logg cvo1 — log<7 spo ), followed by the variation in the mean Fe abundance derived from Fe I lines (column 4), that from 
Fe II lines (column 5), and the average of both (column 6). Similarly, the resulting changes in the K abundances are presented 
in columns 7 (non-LTE K abundance from the 7665 line), 8 (non-LTE K abundance from the 7699 line), and 9 (the average of 
both). Finally, column 10 gives the variation in [K/Fc], which is simply the difference in the results of columns 9 and 6. 



Table 4. Non-LTE reanalysis of Cayrel et al.'s (2004) data. 



Star 


T c s 


logs 


Vt 


[Fe/H] 


EW 7m5 [K/H]ree 


a N LIU 

5^7665 


E W7699 [K/H] 7699 


a N LIE 

^7699 


[K/Fc] 




(K) 


(cm s 1 ) 


(km s" 1 ) 


(dex) 


(mA) 


(dex) 


(dex) 


(mA) 


(dex) 


(dex) 


(dex) 


HD 2796 


4950 


1.5 


2.1 


-2.47 


79.8 


-2.08 


-0.39 


49.9 


-2.13 


-0.30 


+0.36 


HD 122563 


4600 


1.1 


2.0 


-2.82 








44.8 


-2.44 


-0.26 


+0.38 


HD 186478 


4700 


1.3 


2.0 


-2.59 


93.0 


-2.10 


-0.43 


56.0 


-2.22 


-0.29 


+0.43 


BD+17°3248 


5250 


1.4 


1.5 


-2.07 


91.0 


-1.67 


-0.66 


64.9 


-1.71 


-0.46 


+0.38 


BD-18°5550 


4750 


1.4 


1.8 


-3.06 


41.5 


-2.66 


-0.27 


24.0 


-2.68 


-0.24 


+0.39 


CD-38°245 


4800 


1.5 


2.2 


-4.19 
















BS 16467-062 


5200 


2.5 


1.6 


-3.77 


5.6 


-3.35 


-0.20 


2.0 


-3.52 


-0.19 


+0.33 


BS 16477-003 


4900 


1.7 


1.8 


-3.36 








10.3 


-2.98 


-0.21 


+0.38 


BS 17569-049 


4700 


1.2 


1.9 


-2.88 


50.0 


-2.57 


-0.28 


35.0 


-2.51 


-0.26 


+0.34 


CS 22169-035 


4700 


1.2 


2.2 


-3.04 








23.5 


-2.73 


-0.23 


+0.31 


CS 22172-002 


4800 


1.3 


2.2 


-3.86 


5.8 


-3.61 


-0.19 


4.8 


-3.40 


-0.20 


+0.35 


CS 22186-025 


4900 


1.5 


2.0 


-3.00 


39.1 


-2.59 


-0.26 








+0.41 


CS 22189-009 


4900 


1.7 


1.9 


-3.49 


12.5 


-3.18 


-0.20 








+0.31 


CS 22873-055 


4550 


0.7 


2.2 


-2.99 








33.3 


-2.64 


-0.24 


+0.35 


CS 22873-166 


4550 


0.9 


2.1 


-2.97 








37.5 


-2.57 


-0.25 


+0.40 


CS 22878-101 


4800 


1.3 


2.0 


-3.25 
















CS 22885-096 


5050 


2.6 


1.8 


-3.78 


5.0 


-3.51 


-0.19 


2.5 


-3.53 


-0.19 


+0.26 


CS 22891-209 


4700 


1.0 


2.1 


-3.29 


30.2 


-2.88 


-0.23 


16.7 


-2.89 


-0.22 


+0.41 


CS 22892-052 


4850 


1.6 


1.9 


-3.03 








19.8 


-2.70 


-0.23 


+0.33 


CS 22896-154 


5250 


2.7 


1.2 


-2.69 


33.1 


-2.43 


-0.27 


22.2 


-2.37 


-0.25 


+0.29 


CS 22897-008 


4900 


1.7 


2.0 


-3.41 


15.5 


-3.08 


-0.21 


9.6 


-3.02 


-0.21 


+0.36 


CS 22948-066 


5100 


1.8 


2.0 


-3.14 








11.2 


-2.80 


-0.22 


+0.34 


CS 22949-037 


4900 


1.5 


1.8 


-3.97 


2.8 


-3.86 


-0.18 


1.6 


-3.81 


-0.18 


+0.14 


CS 22952-015 


4800 


1.3 


2.1 


-3.43 








10.5 


-3.05 


-0.21 


+0.38 


CS 22953-003 


5100 


2.3 


1.7 


-2.84 








15.6 


-2.65 


-0.23 


+0.19 


CS 22956-050 


4900 


1.7 


1.8 


-3.33 








7.4 


-3.13 


-0.20 


+0.20 


CS 22966-057 


5300 


2.2 


1.4 


-2.62 


34.0 


-2.38 


-0.27 


21.7 


-2.35 


-0.26 


+0.26 


CS 22968-014 


4850 


1.7 


1.9 


-3.56 


6.7 


-3.50 


-0.19 


5.1 


-3.34 


-0.20 


+0.14 


CS 29491-053 


4700 


1.3 


2.0 


-3.04 


48.1 


-2.62 


-0.28 








+0.42 


CS 29502-041 


4800 


1.5 


1.8 


-2.82 


52.9 


-2.46 


-0.29 


36.1 


-2.42 


-0.26 


+0.38 


CS 29502-042 


5100 


2.5 


1.5 


-3.19 


13.9 


-3.00 


-0.21 








+0.19 


CS 29516-024 


4650 


1.2 


1.7 


-3.06 
















CS 29518-051 


5200 


2.6 


1.4 


-2.69 


30.0 


-2.53 


-0.25 


25.0 


-2.34 


-0.25 


+0.26 


CS 30325-094 


4950 


2.0 


1.5 


-3.30 


30.0 


-2.69 


-0.24 








+0.61 


CS 31082-001 


4825 


1.5 


1.8 


-2.91 

















Note. 

The atmospheric parameters (columns 2-5) as well as the EW values for the 7665 and 7699 lines (columns 6 and 9) are from 
Cayrel et al. (2004). Otherwise, the meanings of the given results are the same as in table 1. 



